Spectroscopy of Valley Splitting in a Silicon/Silicon-Germanium Two-Dimensional 

Electron Gas 
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The lifting of the two-fold degeneracy of the conduction valleys in a strained silicon quantum well 
is critical for spin quantum computing. Here, we obtain an accurate measurement of the splitting 
of the valley states in the low-field region of interest, using the microwave spectroscopy technique of 
electron valley resonance (EVR). We compare our results with conventional methods, observing a 
linear magnetic field dependence of the valley splitting, and a strong low-field suppression, consistent 
with recent theory. The resonance linewidth shows a marked enhancement above T ~ 300 mK. 

PACS numbers: 73.21.Fg,78.70.Gq,78.67.De 



The term "valley physics" refers to the study of degen- 
erate valleys in the conduction band of an indirect gap 
semiconductor such as silicon. Valley physics has be- 
come a focal point in the field of silicon spintronics and 
quantum information processing because of the couplings 
between valley and spin states. For instance, in the Kane 
quantum computer |l[ , the interactions between spins are 
strongly modulated by interference between the different 
valleys 2\. Similar concerns exist for spin qubits in a 
quantum well Q. 

Although valley physics has emerged as an important 
field of study, many important experimental questions 
remain unsettled, due to the dearth of valley-sensitive 
measurement techniques, particularly for a two dimen- 
sional electron gas (2DEG). One example is the dis- 
crepency between theory and experiment for the mag- 
nitude of the energy gap between the ground and excited 
valley states (the so-called valley splitting). While the- 
ory predicts that the valley splitting should be of order 
1 meV for a silicon/silicon-germanium quantum well 
experimental measurements can be 10-100 times smaller 
SH, 0, B Q • Typical detection techniques involve beat- 
ing in Shubnikov-de Haas measurements or activation 
energy analyses. These methods are difficult to apply 
with high precision, and they do not work well at the 
low fields of interest for quantum devices. Since valley 
splitting must be large enough to minimize excitations 
outside the qubit Hilbert space of a spin-based quantum 
computer |lOj, it is crucial to understand the low-field 
valley physics, and to perform accurate low-field mea- 
surements. 

In this Letter, we apply the high precision microwave 
resonance techniques developed for spin excitation to 
the problem of valley splitting in a 2DEG. Because of 
the small number of electrons and the low-temperature 
requirements, conventional microwave absorption spec- 
troscopy techniques (which require almost 10 12 electrons 



for an appreciable resonance signal |ll|) are difficult to 
employ in these structures. On the other hand, trans- 
port measurements are naturally suited for probing elec- 
trical characteristics of narrow channels containing rela- 
tively few carriers. Previous studies have combined these 
techniques in the form of electrically detected electron 
spin resonance (ED-ESR), thereby enabling the mea- 
surement of Zeeman splitting in both gallium arsenide 
and silicon 0, GJ 2DEG struc- 
tures. Here, we show that electronic transitions can also 
be driven between the two lowest valley states using mi- 
crowaves to achieve electrically detected electron valley 
resonance (ED-EVR). The advantage of this technique is 
that it allows accurate and dense data acquisition over 
more than a decade of low magnetic fields. 

The Si/SiGe heterostructures used in these experi- 
ments were grown by ultrahigh vacuum chemical vapor 
deposition [l8|. In each case, the 2DEG is located atop 
80 A of strained Si grown on a strain-relaxed Sio.7Geo.3 
buffer layer. The 2DEG is separated from the phosphorus 
donors by 140 A of Sio.7Geo.3. The donors lie in a 140 A 
layer of Sio.7Geo.3 with a 35 A Si cap at the surface. Fur- 
ther details about the structure can be found in reference 
0. Two 2DEGs (51 and 52) were measured at 0.25 K, 
obtaining the electron densities n = 4.2 x 10 11 cm~ 2 
(51) and n = 5.5 x 10 n cm~ 2 (52), and the mobilities 
50,000cm 2 /Vs (51) and 180, 000 cm 2 /Vs (52). 

A schematic of the experimental set-up is shown in 
Fig. [TJ A double lock-in technique is used to measure 
the change in resistance AR XX of the 2DEG as a func- 
tion of the perpendicular field B, in the presence of mi- 
crowaves. Lock-in 1 provides a bias current ranging from 
100 nA to 250 nA, modulated at 701.3 Hz. Lock-in 2 is 
used to modulate the microwave amplitude with 100% 
modulation at 5.7 Hz. The output of Lock-in 1 is fed 
into Lock-in 2, which measures AR XX . Microwaves are 
produced by an HP83650A synthesizer, and are carried 
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FIG. 1: Schematic of the experimental setup, using a double 
lock-in technique (Li and L2). 



sen 

FIG. 2: Valley splitting data for two different samples. Re- 
sults for sample SI (squares) were obtained using EVR. An 
activation energy analysis was used for sample S2 (circles). 
The line represents a linear fit to the transport data for SI. 
Inset: Shubnikov-de Haas data shows splitting of the orbital 
oscillations (O), due to the lifting of spin (S) and valley (V) 
degeneracies. 



down to the sample using a low loss coaxial line termi- 
nating about 5 cm from the surface of the sample in a 
loop antenna. The base of a resonant cavity is replaced 
with a sample stage. The microwave power at the sample 
has a strong frequency dependence because of the open 
cavity and the impedence mismatches along the length of 
the coaxial line. Because of this non-uniformity, a wide 
range of powers (10 /LtW-10 mW) are used to ensure opti- 
mal power delivery. The magnetic field is produced by a 
superconducting magnet and all measurements are car- 
ried out in an Oxford Instruments 3 He cryostat with a 
base temperature of 0.25 K. 

The same experimental set-up can be used to detect 
both ESR and EVR signals. Although we do not re- 
port on ESR here, we observe typical resonances, with 
linewidths on the order of 5 G for SI and 2 G for S2. 
The EVR transition is slightly different than ESR be- 
cause it is not driven by magnetic fields. (The two low- 
lying valley states are orthogonal and unaffected by the 
spin operator, causing the Zeeman transition matrix ele- 
ment to vanish.) However, the electric dipole transition is 
allowed by general symmetry considerations [20j ], which 
also apply to the quantum well geometry. We can esti- 
mate the magnitude of this valley excitation using a one- 
dimensional tight binding (TB) method ,4j. Since the 
two valley states differ only along the z direction (the 
direction perpendicular to the quantum well), only the 
E z component of the microwave field may induce tran- 
sitions. The resulting dipole matrix element is rather 
small. Nonetheless, it is the dominant transition mecha- 
nism. By further positioning the sample at B or E nodes 
in the resonant cavity, it may be possible to drive valley 
or spin excitations selectively, although we do not per- 
form such experiments here. 



Shubnikov-de Haas oscillations provide a rough esti- 
mate of the valley splitting, as demonstrated in the inset 
of Fig. [2] for sample S2 . Both spin and valley features can 
be observed in the data [8|. By increasing the magnetic 
field, we observe a sequential removal of the spin and val- 
ley degeneracies, as indicated by the appearance of split 
spin and valley peaks. In the main figure, we present 
a measurement of the valley splitting (circles), based on 
an activation energy analysis [5j. The analysis provides 
results only at higher magnetic fields. The larger error 
bars reflect uncertainties in the fittingof the activation 
data, similar to previous experiments [5|]. 

In contrast, the electrically-detected EVR measure- 
ment provides narrow error bars, and covers a wide range 
of magnetic fields. Here, we obtain data from 0.27 T to 
3 T on sample SI, with narrow error bars throughout. 
Some typical resonances are shown in Fig. To ana- 
lyze the resonance features, we fit the data. First, the 
background resistance is removed by fitting to a second 
degree polynomial away from the main peak. We find 
that gaussians provide the best representation of the in- 
dividual peaks, with peak widths on the order of 20-25 G. 
Typically, the resonance features account for about one 
part in 10 4 of the total resistance signal. In Fig. [3[ the 
the peak heights have been scaled to unity. 

The fitted peak positions are plotted in Fig. [2] 
(squares), as a function of the perpendicular magnetic 
field. To estimate the error bars, we note that the mi- 
crowave power dependence of the resonance peak shows 
a shift towards higher fields with increasing power. The 
largest observed shift is about 100 G over two orders of 
magnitude in the power. In the figure, we determine our 
error bars using a more generous estimate of 500 G. The 
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FIG. 3: Electrically detected EVR signals at two different 
microwave frequencies with Gaussian fits. 
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resulting data are strikingly linear, with regression giv- 
ing a slope of (24.7 ± 0.4) ^eV/T, and a y-intercept of 
E v (0) = (1.5±0.6) /xeV. A separate linear fit of the valley 
splitting data from 5*2 suggests a slightly larger slope for 
E V (B). Although this fit has larger error bars than the 
EVR analysis, the results appear consistent with theo- 
retical expectations that the valley splitting should scale 
with the 2DEG density as E v oc n [2l| . 

Figure [4] shows the effect of temperature on the res- 
onance peaks. We use the frequency (18.11 GHz) and 
magnetic field (3.0 T) that give the largest valley split- 
ting. Data sets are taken as the sample temperature is 
increased from 0.23 K to 0.35 K and decreased back to 
0.23 K at a slow rate, giving the sample enough time to 
equilibriate. The error bars for Fig. HJb) are obtained 
from the hysteresis in peak widths during the tempera- 
ture cycle. Over this narrow temperature range there is a 
sudden, rapid (seven-fold) increase in the linewidth. Res- 
onances at higher temperatures are difficult to observe, 
although valley splitting can still be observed in the ac- 
tivation energy analysis, for temperatures up to 1 K. In 
addition to thermal broadening of the resonance, we also 
observe a small but reproducible increase in the peak po- 
sition of about 400 G, or 1.5%, indicating a small thermal 
enhancement of the valley splitting. 

While the data in Fig. [2] are internally consistent, they 
give values for valley splitting that are far smaller than 
the theoretical estimates. In the conventional theory of 
valley splitting, the degeneracy of the conduction val- 
leys is broken by the sharp confinement potential of 
the quantum well, which couples the valleys in k space 
USUI- Tight binding and non-equilibrium Green's 
functions techniques obtain estimates of 0.1-1 meV for 
the valley splitting, depending on the quantum well width 
and the internal electric field associated with modula- 
tion doping Q. It has been suggested that the mag- 
netic field dependence arises from the enhancement of 
exchange coupling due to electron-electron interactions 
in the 2DEG " 
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FIG. 4: Temperature dependence of the valley resonance, 
(a) Normalized EVR signals from 230 mK (left) to 350 mK 
(right), (b) Temperature dependence of the corresponding 
fwhm linewidths. 



addition, the linear dependence of E V (B) is inconsistent 
with the expected scaling for the many-body theory [25| . 

Ando has proposed an alternative explanation for the 
suppression of valley splitting 24, 26]. In this picture, the 
sharp confinement potential is still the cause of the valley 
splitting. However, one must also include the effects of 
substrate miscut and rough growth surfaces. Indeed, the 
commercial substrates used for SiGe heterostructures are 
often purposely miscut. The samples used in this work 
were miscut at a 2° angle. Quantum wells grown on 
such substrates will be misaligned with respect to the 
crystallographic axes. For rough surfaces, there will be 
an additional, locally varying misalignment. 

A theory of valley splitting on a stepped quantum well 



22, 24, 2 



However, the many-body ef- 
fect is expected to enhance the valley splitting, in con- 
trast with the suppression observed in experiments. In 



is given in Ref. [2l|, based on effective mass theory. A 
number of the experimental features in Fig.[5]appear con- 
sistent with the theoretical predictions. We now briefly 
discuss the theory, and its implications for our work. The 
valley splitting can be expressed as a simple integral, 
E V =2\J e~ 2tk ° z \F(r)\ 2 V v (r)d 3 r\, where k is the posi- 
tion of the z valley minimum, F(r) is the effective mass 
envelope function, and V v (r) is a valley coupling interac- 
tion, caused by the sharp interface of the quantum well. 
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Each step in the quantum well makes a contribution to 
the integral with the phase e~ 2lk ° Zj , where Zj is the inter- 
face position of the jth step. Since the difference in the 
phase angles on consecutive steps is 2kob ~ 0.85-7T, where 
b = 1.35 A is the atomic step height, the step contribu- 
tions interfere destructively. Thus, an electronic wave- 
function covering many steps will have a valley splitting 
that is strongly suppressed compared to the flat quantum 
wells of Ref. [J]. Since stepped surfaces are ubiquitous 
in conventional semiconductor heterostructures, so too 
is the suppression of the valley splitting. In a magnetic 
field, on the other hand, the electron is confined to a 
finite number of steps, thereby limiting the destructive 
interference. For very large magnetic fields, the electron 
may be confined to a single step. In this limit, the val- 
ley splitting is restored, and E v achieves its theoretical 
upper bound. 

The arguments given in Ref. [2l| for the strong sup- 
pression of E v at very low fields are plausible and con- 
sistent with our data. The partial lifting of this suppres- 
sion due to step disorder and other fluctuations is also 
plausible. In Ref. 2l|, simulations were performed on 
disordered steps, including bunched steps. These obtain 
valley splitting results very close, quantitatively, to the 
data of Fig. [2l However, the shape of the E V {B) curves 
obtained from the simulations depends on the particular 
disorder model, making a definite theoretical comparison 
difficult. Most noteably, the simulation results do not ap- 
pear completely linear, contrary to our experimental ob- 
servations. A particular "plateau" model was suggested 
in that work, as an example of a disorder model produc- 
ing a linear E V (B). However, experimental verification 
of such behavior is not yet available. On the other hand, 
the proposed model of envelope function oscillations at 
B = gives an estimate for the valley splitting which is 
very close to the extrapolated of value for E v (0), obtained 
from Fig. d 

There are other open questions in the valley splitting 
theory of Ref. 2l| . The EVR experiments summarized in 
Figs. 2 and 3 exhibit sharp resonance peaks, indicating a 
well-defined valley splitting. However, models involving 
disorder suggest that electrons can become localized in 
the valley splitting landscape. In this picture, electrons 
may fill both shallow and deep pinning sites, consistent 
with a range of valley splittings. One might therefore ex- 
pect a broad resonance peak, in contrast with experimen- 
tal observations. It is tempting to attribute the marked 
thermal broadening of the EVR peaks to the increased 
occupation of higher energy states. A more complete the- 
ory must also take into account the fact that our electrical 
detection method is inherently dynamical. 

Finally, we discuss the importance of our results for 
quantum computing in a silicon 2DEG. To minimize the 
excitation of the valley states in a spin-based quantum 
computer, the valley splitting should be much larger than 
the temperature. For a system cooled to 100 mK, a 



valley splitting of 100 fieV should be adequate. We 
have noted that a strongly confined electron will ex- 
hibit a larger valley splitting. For the plateau model 
of Ref. [lH, the linear dependence of E V (B) can be ex- 
pressed through E v — C/itR 2 9 2 , where 8 is the miscut 
angle, and R = \/21b = \J1hj\eB\ is the rms radius of 
the magnetically confined wavefunction. We can extend 
this scaling to confinement produced by gate potentials. 
For sample 5*1, with 9 = 2°, our experimental results 
lead to C = 4.1 x 10~ 19 eVm 2 . If we now take R to be 
the quantum dot radius and require E v = lOO^eV, we 
obtain the relation R9 = 36 nm. To attain E v = 100 [ieV 
on a 2° miscut therefore requires a fairly small dot radius 
of 18 nm. However, for a 0.3° miscut, a larger 120 nm 
dot can be used. 

In conclusion, we have performed microwave spec- 
troscopy of the conduction valley states in a strained sil- 
icon 2DEG using transport measurements. Our method 
provides a measurement of the valley splitting over a 
wide range of relatively low magnetic fields (0.3-3 T). We 
also obtain temperature dependent measurements show- 
ing a strong change in the valley splitting behavior above 
300 mK. We compare our data to a recent theory of valley 
splitting on a stepped substrate obtaining quantitative 
agreement, but leaving some open questions. 
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